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We calculate the dc Josephson current Ij for two types of superconductor-ferromagnet (S/F) 
Josephson junctions. The junction of the first type is a S/F/S junction. On the basis of the 
Eilenberger equation, the Josephson current is calculated for an arbitrary impurity concentration. 
If hr -C 1 the expression for the Josephson critical current I c is reduced to that which can be obtained 
from the Usadel equation {h is the exchange energy, r is the momentum relaxation time). In the 
opposite limit lir > 1 the superconducting condensate oscillates with period Vp/h and penetrates 
into the F region over distances of the order of the mean free path I. For this kind of junctions we 
also calculate Ij in the case when the F layer presents a nonhomogeneous (spiral) magnetic structure 
with the period 2n/Q. It is shown that for not too low temperatures, the 7r-state which occurs in 
the case of a homogeneous magnetization (Q = 0) may disappear even at small values of Q. In this 
nonhomogeneous case, the superconducting condensate has a nonzer o triplet component and can 
penetrate into the F layer over a long distance of the order of £t = \J D/2-kT. The junction of the 
second type consists of two S/F bilayers separated by a thin insulating film. It is shown that the 
critical Josephson current I c depends on the relative orientation of the effective exchange field h of 
the bilayers. In the case of an antiparallel orientation, I c increases with increasing h. We establish 
also that in the F film deposited on a superconductor, the Meissner current created by the internal 
magnetic field may be both diamagnetic or paramagnetic. 



I. INTRODUCTION 



Interplay between ferromagnetism and superconductivity in layered structures has attracted a great interest in the 
last years. In a rough approximation, these states are antagonistic to each other and the ferromagnetism, being usually 
much stronger than superconductivity, is supposed to destroy the latter. However, in many cases the coexistence of 
these two phenomena is possible, even if the superconducting critical temperature T c is by order of magnitude lower 
than the Curie temperature of the ferromagnet. Such is the case when dealing with superconductor-ferromagnet (S/F) 
hybrid structures. In these systems the mutual interaction of these two states may lead to significant changes of the 
thermodynamic and transport properties. 

In particular for S/F/S systems in equilibrium, one of the most interesting effects is a phase shift by tt between weakly 
coupled superconductors, the so called 7r-state. T±i£ possibility of the tt -state in S/F/S structures was first predicted 
by Bulaevskii et al.oB and studied in later worksao.The transition to the 7r-state manifests itself in a nonmonotonic 
(and even oscillatory) thickness dependence either of the superconducting critical temperature T c or of the critical 
current J c , and in the change of sign of I c if the exchange field h exceeds a certain value in the S/F/S junction 
(RefsJjl3). Although some experiments on the thickness dependence of T c in S/F structures show that for a certain 
thickness of the F layer, the ground state of the system may correspond to the 7r-phase shift between the adjacent 
superconductors (see e.g. Ref.El), thisikind of coupling was not observed in other experiments(see, for example, Ref.u). 
Only recently, the experiment of Refill on the measurement of the current I C (T) in the Nb/Cu^Nii-^/Nb Josephson 
junction demonstrated unambiguously the transition from the to the 7r-phase difference between the superconductors. 

In all theoretical worksoQ, calculations were performed either in the diffusive limit, in which the Usadel equation 
was applicable, or in the pure ballistic limit, where the elastic scattering by impurities was completely neglected. At 
the same time, very often the parameters characterizing the samples in experiments, such as the sample size, the 
mean free path or the strength of the exchange field, do not correspond to these limits. Therefore, there is a certain 
need to study the Josephson current in the S/F/S structures not only in extreme limits but also in the intermediate 
region of the parameters. 

In this work, we calculate the critical Josephson currents^ in a S/F/S junction for arbitrary impurity concen- 
trations. Since the approach based on the Usadel equationcil (dirty limit) is valid only if the parameter hr is small 
(h is the exchange field of the ferro magn et and r is the momentum relaxation time), we use in an arbitrary case 
the more general Eilenberger equationli§t!j in which, generally speaking, the elastic collision integral is not neglected. 
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As mentioned above, in real experiments the parameter hr may take different values depending on the sample and 
therefore our theory can serve as a good description of the experiments. 

Moreover, in all theoretical works mentioned previously, it was assumed that the magnetic ordering in the ferromag- 
netic layers was homogeneous. However, ferromagnetic materials exhibit generally more complex magnetic structures. 
In strong ferromagnets, like Fe or Ni, the magnetic ground state consists of homogeneously magnetized domains with 
different relative orientations. Also "weak" ferromagnets, like some ternary compounds with a regular lattice of rare 
earths elements, turn out to be superconducting as the crystals undergo a transition into a state with a nonhomoge- 
neous (helicoidal) magnetic order (see Ref.cl and references therein). A similar nonhoipogeneous structure may arise in 
bilayered S/F structures. For example, an experimentliiJ and two theoretical workaiiEj suggested a possible existence 
of a nonhomogeneous magnetic ordering in the ferromagnetic layer in a S/F system. Also in experiments on giant 
magneto-resistance (GMR) in magnetic multilayers employing-,superconducting contacts, nonhomogeneous magnetic 
structures can be created artificially (see e.g. the review Refills and references therein). 

In spite of the importance for the experiments, a theoretical analysis of the influence of a nonhomogeneous magne- 
tization on the properties of S/F junctions is still lacking. Therefore, the second goal of this paper is to investigate the 
influence of nonhomogeneous magnetic configurations on the supercurrent through different kinds of superconductor- 
ferromagnet Josephson junctions. 

In Sec. Ill we consider a S/F/S system, with an nonuniform (spiral) magnetic ordering. We derive an expression for 
the critical current I C (Q) , where Q is the wave vector of the spiral magnetic order. We show that, whereas for Q = 
the transition from the 0-phase state to the 7r-phase state is possible, even small nonzero Q values may restore the 
0-phase state provided that the temperature is not too low. The reason for this is the existence of a triplet component 
of the superconducting condensate in the ferromagnet due to the proximity effect and the nonhomogeneous magnetic 
structure. In the limit hr < 1 this component does not decay over the short distance \jD/h, which corresponds to the 
length of decay of the usual singlet component, surviving up to a much longer distance ~ yjD/2nT (D is the diffusion 
coefficient). The, influence of this triplet component on the transport properties of the S/F mesoscopic structures was 
studied in Ref.113. 

In Sec. IV, we analyze the dc Josephson current in a tunnel junction composed either of two S/F bilayers or of two 
magnetic superconductors. We derive an expression for the critical current I c as a function of the relative angle a 
between the magnetization of both F layers. The most important and surprising result is that for an antiferromagnctic 
configuration, a — n, the current I c increases with increasing exchange field h. The calculated dependence of I c on 
various parameters allows us to make some conclusions not only on the magnetic order of the ferromagnetic materials 
used in S/F,-steuctures but also on nonhomogeneous superconducting states predicted by Fulde, Ferrel, Larkin and 
OvchinnikovNo. 

In Sec.V we show that a Meissner current is induced in the F-region due to the internal magnetic field of the 
ferromagnet. The Meissner current density has a different sign at different points, and the total current in the 
ferromagnet is either diamagnetic or paramagnetic depending on the thickness d of the F film. 

In the Appendix, we present the derivation of the main equations used in this article. All our calculations are based 
on the EilenbergertSE3 or on the UsadeH equations, generalized to the case of a spin-dependent interaction varying 
in space. Another approach based mainly on the Bogolyubov-De Gennes equations was wMely used for the study of 
the spin-injection from a ferromagnet into unconventional superconductors (see, e.g. Ref.E3 and references therein). 
In the present work, we restrict ourselves to the case of conventional superconductors with s-wave pairing. 



II. THE JOSEPHSON CURRENT IN A S/F/S STRUCTURE 

In this section we calculate the dc Josephson current Ij in a S/F/S structure. In order to make the consideration as 
general as possible we use the Eilenberger equationta including the elastic collision term. This allows us to calculate 
I j for an arbitrary impurity concentration and to formulate conditions under which the ballistic or diffusive limits 
can be obtained. In order to find the condensate Green's function f u in the F region in an analytical form, we assume 
that the proximity effect is weak, i.e. <C 1, and linearize the collision term in the Eilenberger equation. This 
assumption can be reasonable for structures with a big mismatch between the Fermi surfaces in F and S, which leads 
to a small transmission coefficient T through the S/F interface. If the coefficient T is of the order of unity, we hope 
that our results are valid at least qualitatively. 

We consider the S/F/S structure shown in Fig. I and assume that the exchange energy h is homogeneous in the 
F region (the case of a non- homogeneous h will be analyzed in the next section). Because of the small interface 
transparency one can neglect the suppression of the order parameter A in the superconductor due to the proximity 
of the ferromagnet. We assume also that there are no spin-flip processes in the ferromagnetic region, i.e. the spin- 
relaxation length is larger than the thickness d of the ferromagnet and there are no spin-processes at the S/F interface. 
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The linearized Eilenberger equation in the Matsubara representation has the following form 

filhdxf + 2 (u m - ih) f = sgnw ( (/) - fj . (1) 

Here T3 is the Pauli matrix, [i = cos#, 8 is the angle between the momentum and the x-axis, / = vpT is the mean 
free path, and u> m — 7tT(2to + 1) is the Matsubara frequency. The angle brackets denote the average over angles: 
< ... >= (1/2) J. d/i(...). Eq. (fj) is complemented by the boundary conditions at x = id/2, which in the case of 
low transparency take the formO 

a =(7/2) sgnwf 3 + s, g s + f s ^ 7Sgnw ( f 3 / s ) , (2) 

L J V / x=±d/2 

where a and s are the antisymmetric and symmetric (with respect to /1) parts of /, 7 = T(/i)/4 is a parameter 
describing the transmittance of the interface, T([i) is the transmission coefficient, g s and f s are the quasiclassical 
normal and anomalous Green's functions of the superconductors. The square brackets denote the commutator. When 
writing the last equality, we neglect the term proportional to /, since |/| ~ 7. The condensate function f s in the 
superconductors can be written as 

f s (±d/2) = (it 2 cos^/2) ± in sin^/2)) f s , (3) 



where f s — A/y/ A 2 + and tp is the phase difference between the superconductors. 

It is convenient to represent the solution of Eq. (|l|) as a sum of the symmetric (s) and the antisymmetric (a) parts, 

/ = § + a (4) 

Substituting this expression into Eq. (|l]) and separating the symmetric and antisymmetric terms, one obtains two 
equations which determine the functions d and s: 

a = -sgnw (pl/K u ) f 3 d x s (5) 

V 2 l 2 dlj - nls = -k u (a) , (6) 

where k u = (1 + 2|w rn |r) — sgnu;2i/iT. Thus, the problem is reduced to finding the solution for Eq. (^J) in the interval 
\x\ < d/2 with the boundary conditions given by Eqs. (Q) and (||). To this end it is convenient to extend formally the 
function s over the whole x— axis and to write Eq. (Jfy in the following form 



2;2o2 - 2 - _ _ 
A* 1 °xx s K u> s ~ K <+ 



(s) + 2^7/ s J2 (^2Cos( ¥ j/2) + (-l) n ir 1 sin(^/2))5(x-(d/2)(2n+l)) 



(7) 



One can prove that the solution of Eq. (^) obeys the boundary conditions. 

Performing the Fourier transformation, we find the solution for the Fourier transform of s 

s k = 2f s .B.F , (8) 

where 

B = (l- K „(M-i))A/ ~ KJ ^ M ~^ - <^ M »] ' ^ 

OO 

F = ^2exp{ikd(2n+l)/2) (if 2 cos(y>/2) + (-l) n in sin(^/2)) (10) 

— OG 

and 

M = ( f dk) 2 + nl . (11) 

The function s determines the dc Josephson current, as well as the variation of the density of states (DOS) due to 
the proximity effect. It is given by the inverse Fourier transformation 

/dk 
— exp(-ikx)s k ■ (12) 
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The current is determined by the expression 



Ij = lG Q N(27n)(2T)TvT 3 £ (a w /x) = (13) 

(jj m 

ic Q iV(2™)(2T)Trf3 ]T ((/i 7 /2) [S(d/2), /.]) , (14) 

where Gq = e 2 /h, N — k%S /tt 2 and 5 is the cross-section area of the junction. In writing Eq. ( |l3] ) we have used the 
boundary condition Eq. (H). The summation over Matsubara frequencies is carried out from m = — oo to m = +oo. 
Substituting Eqs. (|6|-p^|) into Eq. (|l3|), we obtain finally the dc Josephson current 

Ij = I c simp (15) 

where 

I c = G Q NJ C , (16) 

and 



x x f c\k 
J c = 27rTRe ^ ^ f* 2~^ 7M ^ GXP ( ikd ( 2n + ^ 



(17) 



w>0 n=0 



Eq. (|17|) determines the critical current and is valid for any impurity concentration. Its analytical evaluation is 
rather complicated, since it includes summation over Matsubara frequencies, integration over the momentum k and 
the averaging over the angles. Here we will discuss two limiting cases in which Eq. ( p~7[ ) can be simplified. 

a) hr <C L (dirty case). This limit corresponds to a ferromagnet with a weak exchange field h or to an alloy like that 
used in Ref.0, for which the condition h > T c is satisfied. In this case the condition hr <C 1 implies that the quantities 
At and Tt are also small. From Eq. (||) one can determine the coefficient B; in this limit B = [(l 2 /3)(k 2 + fc^)] , 
where n\ — 2{lo — ih)/D and D = w^//3 is the diffusion coefficient. Using Eq. ([Tt]) we obtain for the normalized 
critical current J c 

J c = |( M7 > 2 Re(2^T) £ - , 1 A " 2 . (18) 

4 i! — ' K+i sinh K + a + wf, 

This is the usual expression for the critical current obtained from the Usadel equation in the case of a weak 
proximity effect (cf. Refs.Q0O). In Fig. 2 and Fig. 3 we plot J c as a function of T and d. As has been shown in the 
previous studies, the function J c is a rapidly decaying with T and d function which undergoes several oscillations. This 
oscillatory behavior of the critical current may explain the change from a 0-phase state to a 7r-phase state observed 
in Refill 

We consider now another interesting case. 

b) hr > 1 

This condition corresponds to most of the experiments performed on S/F systems, in which F is a "strong" fer- 
romagnet like Ni or Fe. It does not necessarily mean that we are analyzing the clean case (i.e T c t > 1), since the 
value of the exchange field h can be much larger than T c . For example, if the mean free path I equals ~ 300A, then 
t _1 ~ 300K^> T c , whereas hr > 1 (we take vf = 2.10 7 cm/s). Therefore in the limit ftr > 1 we can deal in principle 
with an arbitrary value of tT c (although realistic materials and samples correspond to the case T c t < 1). It is worth 
mentioning that in this case the use of the Usadel equation is not justified. 

The condition hr ^> 1 implies that k w ^> 1, k u /A1 <§; 1 and, as one can see from Eq. (^|), B = n u {^^,l)/M. 
Performing the integration over k, we find for J c 

J.-^H.E/?( 555CTia? V (1<J) 



uj>Q 



One can see from Eq. ( |l9| ) that the critical current oscillates with varying d or h (cf. Ref.i where the case d <C I 
was considered) and decays with increasing d over the mean free path I. In Fig. 4 and Fig. 5 we plot the dependence of 
J c on d and T calculated numerically from Eq. (|l9|). We see that in this limit the critical current does not oscillate 
with the temperature if the exchange field is temperature independent ( this assumption is quite reasonable in the 
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case of transition metals as Fe or Ni). It can change sign in a hypothetical case T c r > 1. In the limit d/l ^> 1 the 
critical current is exponentially small and one can perform the angle averaging in Eq. (^]). Thus, we find in this limit 

Jc = (ttt) /^ 2 ( 1 ) si £l / T ) ex p (-wx 1 + • ( M ) 

It follows from Eq. (P0[ ) that the critical current oscillates as a function of d or /i and decays with d exponentially if 
d > I and as (hd/vF)^ if d = I. This power law dependence of I c on ft, can be even weaker if 7 depends on /1 (7 has 
a maximum at /1 = 1 and decays with decreasing fx). 



III. S/F/S- JUNCTION WITH AN NONHOMOGENEOUS MAGNETIZATION 

In this section we consider a nonhomogeneous magnetization in the magnetic region. This non-homogeneity can be 
due to domain walls or, as in the case of experiments on giant magnetic resistance (GMR), due to an artificial layered 
magnetic structure. In principle, variation of the magnetic moment on the coordinates can be rather complicated, 
which makes explicit calculations difficult. To simplify the consideration we restrict ourselves with the cases of a 
magnetic spiral structure in the F region with a wave vector Q and calculate the dependence of the critical Josephson 
current I c in a S/F/S junction on the wave vector Q. 

Below we consider only the dirty case (hr < 1) and assume again a weak proximity effect. This means that 7 must be 
small enough: 7 -C \fhr~. In the limit of small hr one can use the Usadel equation for finding the condensate function 
/. However in the case of a rotating magnetization (or equivalently a rotating exchange field h) we need to generalize 
our approach because not only singlet correlators as (V'tV'i) are induced in the F region, but also correlators of the 
type (ipitpi) become nonzero (triplet component). In this case we introduce new 4x4 matrices for the quasiclassical 
Green's functions (see Appendix). The 4x4 condensate Green's function s (to be more exact, its symmetrical part) 
obey the generalized Usadel equation in the Matsubara representation 

-iDg .d 2 xx s+[M h ,s]=0, (21) 

where go = T3 ® do, D = vpl/3 is the diffusion coefficient, and M% = r% ® (aoi\u) m \ + a^hsgnwm cos a) — To ® 
<72hsgnu} n sin a. Here a = Qx is the angle between the z-axis and the direction of h. As in the previous section, we 
neglect corrections to the "normal" Green's function g due to the proximity effect. 

Eq. ( pl| ) is a linear matrix differential equation with space dependent coefficients. This spatial dependence can be 
excluded from the consideration by making a rotation in the Nambu-spin space and by introducing a new matrix s n : 
s = ils n U + , where U — fo <£> <5o cos(a/2) + ir^ ® &\ sin(a/2). After the rotation Eq. ( |2l| ) acquires the form 

d l x Sn - (Q 2 /2) (s„ - As n A + ) + Q (Ad x s n + d x s n A + ) - [f 3 <g> (<7 |u; m | - ihsgnuj m a 3 ) , s n ] + = , (22) 

where A — ir^ <g> d\ and the last term is an anticommutator. If Q = 0, Eq. ( p2j ) coincides with Eq. ( pl| ) and can 
be easily solved. The critical current I c has been calculated in the previous section, and its dependence on T and 
d is presented in Fig. 2 and Fig. 3. It is seen that the critical current I c changes the sign at some h of the order 
of the Thouless energy, td — D/d 2 . The characteristic value of h for the transition to the 7r-state increases with 
increasing T. This result is well known for both types of Josephson (equilibrium and nonequilibrium) junctions in 
which the sign-reversal effect of I c takes place. Earlier than in S/F/S junctions, this effect was observed in 4-terminal 
S/N/S Josephson junctions where a voltage V was-applied between the normal reservoirs and therefore an additional 
dissipative current flows between the N reservoirsc3 were possible. 

The critical current changes sign in these junctions due to a shift of the distribution function in the N electrode with 
respect to the distribution function in the superconductors. In contrast to this case, the sign reversal effect in S/F/S 
junctions is realized at equilibrium conditions. However there is a formal analogy between these two cases because 
the formulae for the critical current, can. be reduced to each other by shifting the energy scale and by replacing h — > 
eV (this analogy was noted in Refs.EaO). 

In the case of a finite Q, the formulae for s become more complicated. In order to make them more transparent, we 
assume that the overlap of the condensate functions s induced by the different superconductors is weak. This means 
that m&x{h, T} should be greater than e<j. Then, we represent the solution of Eq. (p3) for s as a sum of two functions 

§{x) = U [S.s n (d/2 + x).S + + S + .s n {d/2 - x).S] U + . (23) 

The matrix S allows one to take into account the phase difference ip between the superconductors: S — (fa cos(ip/4) + 
i% sin(yj/4))<g> &q. The first and second term in Eq. ( |2^ ) come from the superconductor at x = —(d/2) and x = +(d/2) 
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respectively. The function s n {x) is a solution of Eq. H22 ) for an infinite S/F system with a vanishing phase </? = 0. 
The boundary condition for the new matrix s n has the form 

d x s n + (Q/2) [is„ + s„i+] = -(7/0/. ■ ( 24 ) 

Here f s = if 2 ® 03 / s and 7 = 3(7/i); / s and 7 have been defined in the previous section. It is not difficult to see that 
the function s n (x) has the following form 

s„0) = if 2 ® (ooSoO) + 03S3O)) + in <g> aiT(x) , (25) 

where the functions So(x), Ss(x) and T(x) are the amplitudes of the singlet and triplet component, respectively. All 
these functions may be represented as a sum of three exponentials corresponding to the eigenvalues of Eq. ( |22] ) . For 
example, the expression for S3 (x) is 

S3 (a?) = S34- exp(— k + x) + S3- exp(— K-x) + S31 exp(— mx) . (26) 

Identical formulae may be written for the functions So(x) and T(x) with the factors in front of the exponentials 
denoted as Sq± , Sqi and T± , 7} correspondingly. Analytical expressions for the coefficients and eigenvalues can be 
obtained in the limits of small and large Q, i.e. DQ 2 <C h or DQ 2 h. In the limit of small Q, we find after some 
algebra 



4 = 2(|w m | T sgnuih)/D , nf = 2\u m \/D + Q z (27) 



and 



2 Ms) l( K+K _) Kl ' 53t - (7/s) !(«+«-)«? 



s 3± = tSo± = 7/./2(«±0 .r, = ~(ifs) ^_ : t . ^ = (7/0 . „ 7.2 ^ > ( 28 ) 



where = —2ih/D. 

We note some new important features that appear at finite Q. If Q is zero, only the first two terms in Eq. ( |2^ ) 
are nonzero and the decay is characterized by a short length £p = y/D/h (in case of large enough /i). If Q is finite, 
an additional term (the last term) appears in the formula for S , 3 (x) which, at low temperatures, decays over a much 
larger length of the order y/D/2ivT. Alongside with the last term, the triplet component becomes nonzero. The triplet 
component contains also a long-range term XJ exp(— kix) which increases with increasing Q. Due to the last term in 
Eq. (|26|), the dependence of the critical current on h (or T) is drastically modified even at small Q. 

In order to calculate the current Ij, we can use the general expression Eq. (|l^). In the dirty limit it can be written 
as follows 

I.j = (S/I6p)(2™)(2T)Tr(f 3 CT ) ^scU = (5/I6 ( o)(2^)(2T)Tr (f 3 a ) ( 7 /0 ^ «/«U=d/2 • (29) 

One can show that this expression does not change its form under the [/-transformation and $ in Eq. (|2^) may be 
replaced by s n (d/2). Performing simple operations, we obtain for the current an expression with the same form as 
Eq. @ with 



i c = (s/ P i)f ii 



ui>Q 



exp(~n + d) (Ql) 2 ' 

K+ i + ^F exp( - fi ' d) 



(30) 



The first term in the brackets corresponds to the term [2(re+Z)(sinh K+d)] 1 in Eq. ( p"7j ) in the limit of a large exponent. 
It decays with increasing d over the short characteristic length £p = y/D/h. The second term in Eq. ( |30"| ) is caused 
by the rotation of h along the cc-axis. It decays with d over the characteristic length n^ 1 , which can be much longer 
than £p. Therefore this term leads to a drastic change of the critical current. We calculated numerically the critical 
current I c and presented its temperature dependence for hr — 0.06 in Fig. 6. One can see that for this choice of h 
the critical current is negative at Q = (ir-j unction). However it becomes positive at some finite Ql smaller than hr. 
With increasing Q the long-range term in Eq. (|26| ) and the triplet component increase, reach a maximum and then 
decrease again to zero at large enough Q. 

In the limit of large Q, the coefficients S 3 ±, So± and the triplet components are small. The coefficient S31 — jf s /(kl) 
has the same form as in the absence of h. For the eigenvalues we obtain 
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«± = ±»Q + y/2u m /D , m = ^2(u m /D)+4h 2 /(DQ) 2 . 

We see that in the limit of large Q the solution for / has the same form as in a S/N/S structure (no exchange field), 
i.e. the term (sfs/0) exp(— mx) dominates. The first two terms in the singlet component S 3 (x) which contribute 
to the Josephson current ( see Eq. (pq)) are small. They oscillate rapidly in space and decay over a large distance 
of the order of £t = -JD/2wT (in the limit DQ 2 < h). In the main approximation in the parameter (h/DQ 2 ) the 
temperature dependence of the critical current I c is the same as for a S/N/S junction and we do not present this 
dependence here. 



IV. THE S/F-I-F/S SYSTEM 



In this section we consider a layered system consisting of two F/S bilayers separated by an insulating layer (see 
Fig.l). In this case the Josephson critical current is determined by the transparency of the insulating layer and 
depends on the relative orientation of magnetization in the F layers. 

We assume that the F and the S layers dp,s are thin enough: dp t s < £,f,s, where £p = \J Djh and £s = \J D / A. 
First, we analyze the case of a high S/F interface transparency, i.e. Rs/f < Pf/£,f- Under these conditions all the 
Green's functions are nearly constant in space and continuous across the S/F interface. 

In order to find the Green's functions g R ^ A \ we multiply the components (1,1) and (2,2) of the matrix equation 
( |52| ) (the Usadel equations) by the density of states vp.s in the F and the S layers respectively, and integrate over the 
thickness of the bilayers. Neglecting the influence of one bilaycr on the other (this means that (gd x g) = at the F/I 
interface), we obtain the following equation: 



sgnw 



'03,g 







(31) 



Here the matrix has the same structure as in Eq. (pjl|), but h has been replaced by hp — h (i/pdp) / (vpdp + vsds)\ 
As — A (vsd$) I (ysds + vpdp). We assume that the vector h in the left layer is oriented along the z-axis and has 
the components h (0, sin a, cos a) in the right electrode. One can simplify Eq. ( |3l| ) in the right bilayer with the help 
of the transformation (|5^). In this case one obtains for the both layers the same equation: 



[f 3 (g) (ea + h F a 3 ),g] 



A, 



03, g 



o 



We can solve Eq.(p2|) by making the ansatz 

g = t 3 !g) (a <7 + a 3 <7 3 ) + A s <8> (6o<5o + h^) ■ 

From Eq. ( 

is given by the expression ( |33| ) while in the right bilayer it is given by g^ = U + g V 'U, i.e 

=73(8 (a a + a 3 cos6><7 3 ) - f ® a 3 sin#<7 2 + As <g> (bo cos 9a + b 3 a 3 ) - f 3 A$ ib Q sm9ai 



(32) 



(33) 



and the normalization condition fl5q) one can obtain the coefficients a's and 6's. In the left bilayer g 



According to Eq. (57) only the coefficients bo and 6 3 will enter in the expression for the Josephson current, and they 
are given by 



( b 3)l,r 



1/1 



2 U+ t- 



and (b )i tr 



l,r 



-(- — - 



where £± = y£± — I Asp, and e± = iuj m ±h. By writing As = |As| exp(itp) in the right side one obtains the following 
expression for the critical current 



3V c (a) ee eI c R b = 2nTAiA r ^| Re (r-) Re (7~ 

m>0 ^ V4m/; \4m 



cos Mm ( — J Im( — 



(34) 



where = y (uo m + ihp) + A| and Rb is the tunnel resistance of the I layer. Formula (|34| ) coincides with the 

formula presented in Ref.E3. In the latter work the authors considered a Josephson junction consisting of two magnetic 
superconductors with an oscillating magnetic order. Thus, we have shown that the system of Fig. 7 and that of Ref.c3 
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are equivalent. However the authors of RefBa did not-consider some interesting properties of such structures. We 
note that the same structure was also analyzed in Ref.c3, where the critical current was calculated for different S/F 
interface transparencies. The authors have found the conditions under which the system undergoes a transition to 
the 7T state; however they analyzed only the case of parallel magnetization. 

Here we consider two limiting cases: a) a parallel relative orientation of the magnetizations, i.e. a — and b) an 
antiparallel orientation: a = n. 

In the case a = according to Eq. ([34J), the critical current is given by the expression 

e y cTT = eI cn R b = 4.TAI £ - 2 ^ + %f A 2h2 ■ (35) 

In writing Eq. ( |35"| ) we assumed that hp and | Ag| are the same in both-bilayers (symmetric structure). The dependence 
of the critical current on the exchange field Hf was presented in RefO. At T = the current I c is constant up to the 
value hp = Ao where it drops to zero; Ao is the effective energy gap As at zero temperature and zero exchange field. 
This is a consequence of the fact that the order parameter A is also constant. We do not consider here .a-possible 
transition to the LOFF phase predicted by Larkin and Ovchinnikov (LO)tZl and Fulde and Ferrell (FF)ll3 for the 
region 0.755 Ago < hp. We argue that since the homogeneous superconducting state in this region is a metastable 
state, its realization is possible. Nevertheless our result is definitely valid for the region of small hp, and a possible 
transition to the LOFF would manifest itself in a drop of the the critical current. 

More interesting is the case when the relative orientation of the magnetizations is antiparallel, i.e. a = ir. Then, 
the critical current is given by the expression 

eV cU (n) = eI cU R b = AttTA 2 s V 1 . (36) 

In this case the dependence of I c on hp is completely different from that given by Eq. (|35|). The critical current 
determined by Eq. (|3^) increases with increasing hp (i.e. with increasing either h or dp) and even diverges at zero 
temperature when hp — > Ag. Of course, there is no real divergence of I c since-for example, finite temperatures smear 
out this divergency. The dependence of eV c /Ao on h was presented in Ref.tll. The critical current has a maximum 
at some value oi hp close to Ao. With decreasing T the maximum value of I c increases and its position is shifted 
towards Ao. For arbitrary relative orientations of magnetizations the expression for V c (a) can be presented in the 
form 

V c (a) = V cn cos 2 (a/2) + V cU sin 2 (a/2) . (37) 

Therefore, the singular part is always present and its contribution reaches 100% at a — it. All the conclusions given 
above remain valid also for two magnetic superconductors with uniformhz-joriented magnetizations in each layer. We 
note that in contrast to the case of the spiral structure analyzed in Ref.Ej, no 7r-state appears in our model for any 
effective exchange field hp < Aq (at larger hp the superconductivity is destroyed). As in the previous case of parallel 
orientations, the state with hp — A might be unreachable for the antiparallel orientation due to the appearance 
of the inhomogeneous LOFF state. However the singular behavior of I c can be realized at smaller values of h in 
a structure with large enough S/F interface resistance Rs/f- In this case the bulk properties of the S film are not 
changed by the proximity of the F film (to be more precise the condition Hs 'f > (vpdp/vsds)pF£F must be satisfied; 
pp is the specific resistance of the F film). Then, as one can readily showEll, a subgap e sg — (Dp) F / arises 
in the F layer. The Green's functions in the F layer have the same form as in Eq. ([33]) with As replaced by e sg . The 
singularity in I c (hp) occurs at hp equals to e sg , and the LOFF state does not arise because the subgap e sg is not 
determined by the self-consistency equation. 

For completeness we note that the effect of the relative opesttation of magnetization in the F films on the critical 
temperature T c of the superconductor was analyzed in RefsEirE3 for a F/S/F structure. 

V. ORBITAL EFFECTS 

In the preceding sections we have presented the formulae for the condensate function s (or s) in the ferromagnetic 
regions induced by the proximity effect. The amplitude of s is determined by the interface transparency, i.e. by the 
parameter 7, and the penetration length depends essentially on the parameter hr. The internal magnetic field B 
of the ferromagnet, induces screening currents and leads to some suppression of the condensate function. We have 
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neglected this suppression due to Meissner currents (orbital effects). In order to understand why this approximation 
can be justified we estimate now the magnitude of these effects. 

In the dirty limit the depairing rate due to Meissner current is determined by the energy Dp 2 , where p s = A (d) / 4>o 
is the condensate momentum, A (x) = Bx is the vector potential and 4>q is the magnetic flux quantum. The depairing 
factor can be neglected in the Usadel equation provided the condition 



Dp 2 s < h 



(38) 



is satisfied. For example, for B = lkG, d = 100A, vf — 2.10 7 cm/s and l~rfwe obtain Dp 2 s ~ 50mK. If the condition 
Eq. (|38| ) is met, the condensate function s{x) (to be more exact, its Fourier transform) in the ferromagnet is given 
by Eqs. §0). 

Due to the condensate penetration and the intrinsic magnetic field of the ferromagnet, the Meissner currents arise 
in the F region. In order to analyze this issue in more detail, we consider the S/F system of Fig. 8. 

If we consider the diffusive regime, the condensate function can be found as it was done in section II or directly 
from the Usadel equation. In this limit (hr <§; 1) it has the form: 



s w {x) = 3(/i7)/ s 
The current density is expressed in terms of s(x) as 



cosh K+(rf — x) 
{n+l) sinh(K + <i) 



j(x) = ~<7{Bx/M{2ttT)Rc 



(x). 



(39) 



(40) 



In Fig. 9 we plot the spatial dependence of the current density j(x) which is spontaneously induced in the ferro- 
magnetic film. We can see that j(x) changes sign with varying x. According to the results of section II, in the case 
hr > 1, the current density changes sign many times on the mean free path. 

Integrating the current density j(x) given by Eq. (f40|), we find the total current I = J Q dxj(x) 



I=-\ (3(M7)) 2 ( r(Sd 2 /0o)(27rT)Re X ' 



fs 



^ (n + l) 2 sinh 2 (nj) 



sinh 2 (K+d) 
(k+d) 2 



(41) 



In Fig. 10 and Fig. 11 we presented the total current / as a function of d and hr respectively. It is seen that the total 
current also changes sign, i.e. in the ferromagnetic film either a diamagnetic or paramagnetic current is induced 
depending on the relation between d and 

In the analysis presented here it was assumed that the exchange field h is homogeneous. In a multidomain ferro- 
magnet one expects a more complicated spatial distribution of the Meissner current. 



VI. CONCLUSION 



We analyzed specific features of a supercurrent in superconductor/ferromagnet structures. In section II we have 
calculated the Josephson current Ij in a S/F/S junction. It turns out that the product hr of the splitting energy h 
and the momentum relaxation time r is an important parameter, which determines the approach to bjej^ed in the 
problem. In the dirty limit, i.e. when hr <C 1, Ij can be obtained from the Usadel equation (Refs.QaO). In this 
limit, the change of sign of the critical current I c occurs if the thickness of the F layer d is of the order of £p = \/ ' D jh. 
The condensate function in the F layer decays exponentially over the length £p and undergoes oscillations with the 
same period. In the opposite limit (hr ^ 1) the condensate function oscillates in space with the period vp/h (as in 
the pure ballistic case considered in Reffl) and decays exponentially on the mean free path /. The critical current I c 
decreases with feasa power-law function and is not exponentially small if d ~ I. 

We have also studied the influence of different inhomogeneous magnetic structures on the critical current through 
S/F structures. In section III we considered a S/F/S sandwich with an inhomogeneous magnetic order in the F layer 
described by a vector Q. In the case Q = we obtained the well-known transition from the 0-phase to the 7r-phase 
state. We have also shown that even for small values of Q and not too low temperatures this transition may not take 
place (Fig. 6). The reason for a qualitative change of the I c (hp) is a long-range term in the singlet component of the 
condensate function /. This term arises together with the triplet component if Q differs from zero. The long-range 
part of / decays in the F film on a length of the order WD/2tfT, which can be much longer than the characteristic 
length (~ \J D jh) of the decay of / in a homogeneous F layer (Q = 0). 
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Our results may be applied to ferromagnets containing domain walls and magnetic multilayers with nonhomogeneous 
magnetic structures. We used the quasiclassical Green's function approach to describe such structures in a quantitative 
way. In section IV it was shown that for an antiferromagnetic configuration in the S/F-I-S/F junction, the dependence 
I c (h) shows an anomalous behavior: the critical current increases with increasing h or dp. This means that the 
Josephson critical current in a junction formed by two ferromagnet-superconductor bilayers may be even larger than 
the critical current in a similar Josephson junction S/I/S. In the last section we have considered a S/F bilayer structure. 
The Meissner currents which are spontaneously excited due to the internal field of the ferromagnetic film have been 
calculated. The current density oscillates along the x-axis and the total Meissner current in the F film may be either 
diamagnetic or paramagnetic depending on the thickness d and on the exchange field h. All the effects analyzed in 
our work can be verified experimentally. 
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APPENDIX 

In this part we present some general formulae for superconductivity in the presence of an exchange field. We 
consider structures in which the superconducting pairing and the exchange interaction of electrons with ordered, 
localized magnetic moments take place. The Hamiltonian describing the system under consideration has the form: 

= H { a i 5 pp' + eV ) + ^) 5 "' - (h.(j)] a sV - (Aa+ a+ p , + c.c.) } (42) 

{P,s} 

Here the summation is carried out over all momenta (j>,p') and spins (s, s'), £ p — p 2 /2m — ep is the kinetic energy 
counted from the Fermi energy ep, V is a smoothly varying electric potential, Ui mp = U (p—p 1 ) is a potential describing 
the interaction of electrons with non-magnetic impurities, h is an effective "magnetic field" caused by the exchange 
interaction of spins of the free electrons with spins of the localized magnetic moments. The notation s, p means 
inversion of both spin and momentum. The order parameter A must be determined self-consistently. In order to 
define the Green's function in a customary way we introduce new operators c+ s and c ns , which are related to the 
creation and anhilation operators o+ and a s by the relation (we drop the index p related to the momentum) 



a 8 ,n = 1 
ai,n = 2 . 



The index n operates in the particle-hole (Nambu) space, while the index s operates in the spin space. The operators 
c ns obey the commutation relations 

CnsC n i s t ~t~ C n / s rCns $nn'$ss' i (43) 
CnsC-n' s' ~t~ c n' s' c ns ^rin'&ss 1 • (44) 

In terms of the c ns operators the Hamiltonian can be written in the form 

H = c+ s W( n „/)( ss /)C„/ s ' , (45) 

{p,n,s} 

where the summation is performed over all momenta, Nambu and spin indices. The matrix Ti is given by 

^ = \ { [(£p<W + eV ) + Ui,np] t 3 <g> <t + A <g> <t 3 - h [(f ® <t 3 ) cos a + (f 3 ® <t 2 ) sin a] j (46) 

The matrices fj and cr, are the Pauli matrices in the Nambu and spin space respectively; i — 0,1,2,3, where f 
and (T are the corresponding unit matrices. We have assumed that the exchange field h has the components h = 

h(0, sin a, cos a); this is the case we consider in the next sections. The matrix order parameter equals A = TiReA — 
T2-ZmA. Now we can define the matrix Green's functions (in the Nambu£g>spin space) in the Kcldysh representation 
in a standard way 
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G(t i ,t' k ) = 7 (Tc(c ns (t i )c+A^))), 



(47) 



where the temporal indices take the values 1 and 2, which correspond to the upper and lower branch of the contour 
C, ruauuig from — oo to +oo and back to -co. The quasiclassical Green's functions g{U,t' k ) are defined as usual ( 
RefsJla " 



<?(Pf, r) = - (f 3 (g) a ) / d£ p G(ti, t' k ; p, r) . 



(48) 

We also introduce, as it was done by Larkin and Ovchinnikov (RefS) , a hypermatrix g . The matrix elements of g 
are the retarded g R , advanced g A and the Keldysh (/^component. Thus, g has the form 



g A 



(49) 



The functions g R { A > and g K can be expressed in terms of the time-ordered Green's functions g(ti,tk) as follows 

g R(A) = g(tut[) - 9(ti(2),t' m ) (50) 

g K =g(t u W)+g(t 2 ,t[). (51) 

The equation for the hypermatrix g can be-easily derived in the same way as it was done for the case of a superconductor 
with spin-independent interactions ( RefE9). We are interested in the diffusive limit. The symmetric component of 
the matrix g with respect to the momentum direction pf satisfies the equation: 



A®<7 3 ,g 



- iDV (gVg) + i (f 3 ® a .dtS + d t >g.T 3 ® a ) + eV{t)g - geV(t') 
Here D — vl/3 is the diffusion coefficient, Mh = h (t^ ® 173 cos a + fo (E> 62 sin a) , and 

Eq.The (|5^ ) is supplemented by the normalization condition 

g-g = 1 , 

and by the self-consistency equation 



[M h ,e] =0 



(52) 



A = A Tr (fx - it 2 ) ®a 3 J deg K 



(53) 



(54) 



If the magnetization, and hence the exchange field h, is constant in the ferromagnetic layers, the angle a in Eq. (|5 
can be excluded with the help of the following unitary transformation 



g = U+.g.U 



(55) 



where U = fo^iJo cos(a/2)+i sin(a/2)f3<g)<7i. We consider F/S staictures, therefore we need the boundary conditions 
for Eq. (p2) at the interface between the conductors "a" and "fr"E2l 



ga^a;ga 



Pa 



2R 



a lb 



[Sa 1 ghj j 



(56) 



where p a is the specific resistivity of the conductor "a" and R a /b is the interface resistance per unit area. We assumed 
that there are no spin-flip processes at the interface. In the presence of spin-processes at the boundary the condition 
(|56|) can be generalized (Ref.o). The current density is determined by the usual expression 



I J = ^-Tr (f 3 ® a ) J de (g R .d x g K + g K .8 x g J 



(57) 
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FIG.l: The S/F/S system. 

FIG. 2: Dependence of the normalized critical current J c /j 2 Ao on the thickness d of the ferromagnet for Aor = 0.05 
and T/Ao = 0.01. The curves for hr = 0.1 and hr = 0.5 are multiplied by a factor 10 for clarity. 

FIG. 3: Dependence of the normalized critical current J c /j 2 A n on the temperature for A t = 0.03, hr = 0.06 and 
d/l = IT. 

FIG. 4: Dependence of the normalized critical current J c /j 2 Ao on the thickness d of the ferromagnet for Aqt = 0.05 
and T/A = 0.1. 

FIG. 5: Dependence of the normalized critical current J c /-y 2 A on the temperature for d/l = 1. The two upper 
curves correspond to the case hr = 5, and the the two lower curves to the case hr = 10. For clarity, in the three lower 
curves, the values of the critical current have been multiplied by a factor 5. 

FIG. 6: The dependence of the critical current on T for hr = 0.06, A t = 0.03, d/l = n and different values of Ql. 

FIG.7: SF/I/SF system. 

FIG.8: S/F bilayer. 

FIG.9: Spatial dependence of the current density. Here J = r^-^f, d/l = 10, T/A = 0.1 and A r = 0.05. For 
clarity, the values of the current density for hr = 0.1 have been multiplied by a factor 5. 

FIG. 10: Dependence of the total current on the thickness d of the ferromagnet film for Aor = 0.05 and T/ Aq = 0.1. 
Here 7=^-,^ 

FIG. 11: Dependence of the total current on the parameter hr for A r = 0.05 and T/A = 0.1. Here / = 7 ^ 0ff ^k. 
At hr = 0: 7(0) = -34, -25, -18 for d/l = 1, 5, 10 respectively 
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